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A MATHEMATICAL MODEL FOR THE CUSTOMER 
DYNAMICS BASED ON MARKETING POLICY 

CESAR M. SILVA, SILVERIO ROSA, HELENA ALVES, AND PEDRO G. CARVALHO 


Abstract. We consider a compartmental model to study the evolution of the 
number of regular customers and referral customers in some corporation. Tran¬ 
sitions between compartments are modeled by parameters depending on the 
social network and the marketing policy of the corporation. We obtain some 
results on the asymptotic number of regular customers and referral customers 
in several particular scenarios. Additionally we present some simulation that 
illustrates the behavior of the model and discuss its applicability. 


1. Introduction 

In marketing literature it has been successively referred the importance of cal¬ 
culating the value of a customer. In fact, such indicative value enables firms to 
select those customers that can add profit and consequently constitutes an impor¬ 
tant information to segment the market and efficiently allocate marketing policy 
resources dung. 

The objective of this work is to establish and study a compartmental model, 
mathematically translated into a system of ordinary differential equations, for the 
evolution of the number of customers of some firm, assuming that the customers 
are divided in two subgroups corresponding to different profitabilities. 

Until recently, the value of a customer for a company was based on the present 
value of future profits generated by a customer over the full course of their dealings 
with a particular company, this is the customer life-time-value (CLV) [9]. How¬ 
ever, other authors refer to the importance of including not only the present and 
future revenue from the customer purchases, but also the value of the potential 
to influence other customers under incentives on behalf of the company (customer 
referral value) or by own initiative (customer influencer value) Customer in¬ 
fluencing behaviors consists of the intrinsic behaviors motivating the customer to 
persuade and influence other customers without there being any type of reward on 
behalf of the company and thus designated the customer influencer value (CIV). 
In turn, the patterns of customer recommendation are related to the acquisition 
of new customers due to company initiatives that reward recommendations made 
to other customers, and thereby establishing the customer referral value (CRV). 
According to Kumar et al. [9], these components are mutually interwoven. Thus, 
CLV positively correlates with CRV (although only up to a certain point and in 
an inverted U-shaped relational curve, which means customers reporting average 
CLV are those most interested in company referral programs) and CLV is positively 
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related with CIV (with an inverted U-shaped relationship in effect between these 
two concepts). Much of the literature has focused on the customer referral value 
through the influence customers might have on the formation of other customers 
attitudes (0) in the purchasing making decision {2J and in the reduction of other 
customers perceived risks [5], but little is known about how this processes occur. 
Since customer referral value and influence!' value might have a great impact for 
companies, these latter try to identify the most influential customers [6]. 

A number of studies allow us to think that the customers of a firm can be 
classified into several groups according to their influential role over other potential 
buyers. In imperfect competitive markets information is not purely transparent; 
some persons are more able than others of influencing people to become a customer 
of that firm. It is also acceptable to assume that knowing the referrals among 
each firms customers and quantifying their influence constitutes an important asset 
for the firm competitive advantage, although all customers are important, referrals 
would be more valuable. 

We mostly agree with Marti and Zenou [12j when they state that physics/applied 
mathematics are capable of reproducing many real networks but never reach to 
explain why they emerge; the economists are very precise to explain why they 
emerge but their approach does a poor job in matching real world networks. That 
is why some game theorists are now improving models which take networks as given 
entities and study the impact of their structure on individuals outcomes. 

Based on the network theory some models have been tested to study the way in¬ 
fluential customers can influence other consumers. For instance Kiss and Bichler [6] 
tested real network models, simulated networks and diffusion models to predict in¬ 
fluence between customers based on their position within the network. However, 
as the authors mention this analysis not always is possible if we do not know or 
do not have information regarding the customer social network. Therefore, other 
models are needed to try to explain these processes. 

In this work we propose a model suitable to describe the dynamics of the num¬ 
ber of customers of a given firm. This model is given by a system of ordinary 
differential equations whose variables correspond to groups of customers and po¬ 
tential customers divided according to their profile and whose parameters reflect 
the structure of the underlying social network and the marketing policy of the firm. 
We intend to understand the flows between these groups and its consequences on 
the raise of customers of the firm. We also want to highlight the usefulness of these 
models in helping firms deciding their marketing policy. 

Specifically, the main objectives of our study is threefold: we intend to obtain 
theoretical results concerning the long term behavior of the number of customers 
in various scenarios, we want to present some simulation aimed at illustrating the 
possibilities of application of our model and, finally, we want to discuss the benefits 
and limitations of this type of analysis. 

As referred, we will consider a compartmental model. As far as we are aware, 
this is the first time such type of mathematical model is considered in the context 
of marketing research. We believe that this type of model can be fruitfully explored 
in this context. This believe is based on the fact that compartmental models have 
proved to be an important tool not only in the natural sciences, particularly in 
mathematical epidemiology [3] and in population biology nsma, but also, with 
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increasing notoriety in recent years, in the context of economy and other social 
sciences [EHIHl HUE]. 

We consider a continuous compartmental model with four compartments, rep¬ 
resented by the graph in Figure [l] and governed by an autonomous system of four 
ordinary differential equations. 
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Figure 1. The compartmental model 

Measuring time in years, we consider the following (pairwise disjoint) compart¬ 
ments: R(t), the referral customers in time t, C(t), the regular customers in time 
t, P#(t), the potential referral customers in time t and Pc(t ), the potential regular 
customers in time t. To model transitions between compartments we consider the 
following parameters: Ai, the natural transition rate between Pc and C, given by 
the number of potential regular customers that become regular customers without 
external influence per year over the number of potential customers (by “without 
external influence” we mean without being influenced by marketing campaigns or 
referral customers); A 2 , the referral pull effect, given by the average number of cus¬ 
tomers that a single referral brings (with no additional incentive) per year over the 
number of potential customers; A 3 , the natural Pr to R transition rate, correspond¬ 
ing to the number of potential referral customers that become referral customers 
without external influence per year over the number of potential referral customers; 
m(t), the undifferentiated marketing costs, corresponding to marketing costs asso¬ 
ciated to undifferentiated marketing campaigns per year; A 4 , the pull effect due to 
undifferentiated marketing, corresponding to the quotient of the outcome of undif¬ 
ferentiated marketing campaigns per year by the number of potential customers (by 
“outcome of undifferentiated marketing campaigns” it is meant the number of po¬ 
tential customers that become customers in the sequence of undifferentiated market¬ 
ing campaigns per unitary marketing cost per year); m#(t), the referral associated 
marketing costs, corresponding to marketing costs associated to referral directed 
marketing campaigns per year; A 6 , the pull effect due to referral directed market¬ 
ing, given by the referral directed marketing campaigns outcome over the number 
of potential customers (by “referral directed marketing campaigns outcome” it is 
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meant the average number of additional customers that a single referral can bring 
with incentives per unitary marketing cost per year); A5, the non-central/central 
transition in the social network equal to the number of individuals non-central in 
the social network that become central over the total number of individuals in the 
social network; A7, the central/non-central transition in the social network, given 
by the number of individuals central in the social network that become non-central 
over the total number of individuals in the social network; /?i, the regular customer 
defection rate, equal to the number of regular customers that cease to be customers 
over the number of regular customers; 02, the referral defection rate, given by the 
number of referrals that cease to be customers over the number of referrals; e, the 
customer & potential customer defection rate, corresponding to the number of in¬ 
dividuals that leave the universe of customers and potential customers per year 
over the number of customers and potential customers (by “number of individuals 
that leave the universe of customers and potential customers per year” it is meant 
the number of customers and potential customers that cease to be in the set of 
customers or potential customers per year due to emigration, death, etc.); 7, the 
customer and potential customer recruitment rate, given by the number of individ¬ 
uals that enter the universe of customers and potential customers per year over the 
number of customers and potential customers (by “number of individuals that leave 
the universe of customers and potential customers per year” it is meant the number 
of customers and potential customers that cease to be in the set of customers or 
potential customers per year due to immigration, etc.); a, the referral recruitment 
rate, equal to the number of referrals that enter the universe of customers and 
potential customers per year over number of individuals that enter the universe of 
customers and potential customers per year. 

Our model can be translated into the following system of differential equa¬ 
tions m to be studied along this paper. 


( C' — A7 R — (e + d 1 + A5 )C + (Ai + mX/CjPc -f A2 RPc 
j R' — A5 C — (e + /3 2 + A 7 )R + (A3 + m\i)Pn + (A2 -t- rirXq)RPr 
Pc = (1 — o)l + PiC + A 7 Pr — (e + As + Ai + m\A)Pc — A2 RPc 
, P11 — Oi'y 02 R -f- As Pc — (e 4 - A7 + A3 -t- 7 ti\a)Pr — (A2 + ttirXq^RPr 


Notice that A2 RPr correspond to the average number of referrals that are brought 
with no incentive by the referrals per year, that A2 RPc is the average number 
of regular customers that are brought with no incentive by the referrals per year, 
that mn(t)XeRPR is the average number of additional customers brought due to 
incentives per year, that m(t)\4Pn is the number of potential referrals that become 
referrals in the sequence of undifferentiated marketing campaigns per year and 
that m(t)X4Pc is the number of potential regular customers that become regular 
customers in the sequence of undifferentiated marketing campaigns per year. 

This paper is divided in the following way: in section[ 2 ]we state our main results 
concerning the asymptotic behavior of the number of regular customers and referral 
customers, in section[ 3 ]we present some simulation with the objective of illustrating 
our theoretical results, in section[ 4 ]we prove our results and, finally, in section[ 5 ]we 
discuss the results obtained. 
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2. Main Results 

One of the first natural issues to address when studying a compartmental model 
is the existence and stability of equilibrium solutions. We obtain several results on 
the existence and stability of equilibrium solutions in of model © in this section. 
We first derive an auxiliary result. Given <5^0 define the sets 

Aj = {(a;, y, z , w) G (R + 4 : y/e - <5 < a; + y + 2 + w < j/e + <5}. 

We have the following result that shows that the total population in system © 
converges to the ratio 7 /e, independently of the nonnegative initial conditions con¬ 
sidered. 


Lemma 1 (Asymptotic behavior of the total population). Let e > 0 and let 
(C(t ), R(t ), Pc(t ), Pn(t)) be some solution of system © with nonnegative initial 
conditions: C(to) ^ 0, R(t 0 ) > 0, P(t Q ) > 0, P(to) ^ 0. Then: 

a) for all t ^ to, we have C(t), R(t), Pc(t), Pn(t) ^ 0; 

b) we have lim C(t) + R(t) + Pc it) + PR.it ) = —. In particular, given 5 > 0 

t —>+00 £ 

and any solution {C(t), R(t), Pc(t), Pn(t)), with nonnegative initial conditions 
there is T > 0 such that (C(t), R(t), Pc(t), Pfi(t)) G As for all t ^ T and any 
equilibrium solution is in the set Ao = {(x, y, z, w) G (Rq ") 4 ■' x+y+z+w = 7 / 5 }. 

Remark 1. The case e = 0 is not a very interesting case since it corresponds to 
the situation where there is no customer & potential customer defection rate which 
is not a realistic assumption. Nevertheless it is easy to check that, if £ = 0, then, 
given initial conditions Co, Rq, Pep an d Pr.o, we have 

lim C(t) + R(t) + Pc(t) + Pnft) = 7 1 + Co + Ro + Pc,o + Prp- 

£—>•+00 

In particular if e = 7 = 0, then the total population remains constant. 


We now obtain a result on the existence of equilibrium solutions. Under the 
assumption of positivity of the defection rate, the referral pull effect and the non 
central/central transition in the social network we conclude that there are one, two 
or three equilibrium solutions, depending on the number of real roots of some third 
degree polynomial. We need to define the constants 


( 2 ) 

and also 


_ 7(Q£ + As) 

e(e + A 5 + A 7 ) 


and 


7((1 ~ Q+ + A 7 ) 
e(e + A 5 + A 7 ) 


u — £ ~\~ /?i A 5 Ai A 171X4 and v — £ ~\~ /?2 d - A 7 A 3 A 7x1X4. 


Theorem 1 (Equilibrium solutions). Let e, A 5 > 0. Then: 

a) system © has up to three equilibrium solutions (R *, C *, P^, Pff). The first 
component, R*, is always a nonnegative solution of the cubic equation 

(3) aR 3 + bR 2 + cR + d = 0, 

where a = —X 2 (X 2 + tur A 6 ), b = A 2 (A 2 + to_rA 6 )p — u(X 2 + mRX 6 ) - X 2 v, 
c = X 2 (X 3 +mX 4 )p + u(X 2 + mRX 6 )p+{X 7 + X 2 q)X 5 — uv and d = (X 3 +mX 4 )pu + 
(Ai + mA 4 )gA 5 ; 
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b) any equilibrium solution (R *, C *, P R , Pq) is obtained in the following way: R* 
is a nonnegative solution of 0 and PJ) = p — R*, P £ = q — C* and 

n* — ( e + $2 + A t)R* — (A3 + toA 4 + (A2 + m R \o)R*){p — R*) 

(, 4 J 0 —-1-) 

A5 

are nonnegative constants. 

In the next result we discuss the asymptotic behavior of solutions of 0 under 
some assumptions on the parameters that roughly correspond to require that the 
referral pull effects are bounded by some functions that we can identify with the 
other “forces” in the model such as the natural transition rates, the pull effects 
due to undifferentiated marketing and the defection rates (see equation 0 ). In the 
following theorem we were able to show that, under the mentioned assumptions, the 
asymptotic behavior of the solutions of 0 can be obtained by the two dimensional 
autonomous system (Ell- 


Theorem 2 (Asymptotic behavior of solutions). Let e > 0 and assume that 
(5) 

f 2 e -f- 2 /?i + A5 + 2 Ai + 2 mA 4 2 s -t- 2/I2 T A7 + 2A3 + 2 mA 4 + 2(A2 m R X§)p 1 

min < ---7-, ---r- > > 1. 

( A7 + q\2 A5 + q\2 J 

Consider the system 

,g, j C'a = XjRa — (e + / 3 l + A5)Ca + (Al + IT1A4 + \2Ra)(q — Co) 

( R'a = X§C a — (e + P2 + A7 )Ra + (A3 + mA 4 + (A2 + mnXg)R a ) (p — R a ) 

and set Pr , a (t) =p — R a (t) and Pc, a (t) = Q — C a (t). Then the asymptotic behavior 
of C, R, Pc and Pr in system 0 is the same as the asymptotic behavior ofC a , R a , 
Pc,a and P Rta . Namely if (C(t), R(t), Pc{t), Pn{t)) is a solution of 0 with initial 
condition {C(t 0 ), R{t 0 ), P c {t 0 ), P R {t 0 )) = {C 0l R 0 , Pc,o, Pr,o) and (C a (t), R a (t)) is 
a solution of with initial condition (C a (to), R a (to)) = (Cq,Rq) then 


lim | C(t) - C a {t) | + |I?(f) - R a (t)| + I Pc{t) - Pc,ait) | + | P R (t) - P R ,a(t)\ = 0. 

£—>■+00 

In the next two theorems we discuss two particular situations where we analyse 
the existence of equilibriums and their stability. The vector fields plotted with the 
objective of illustrating the situations correspond to the reduced system 0 but we 
considered situations where Theorem [ 2 ] apply so that the asymptotic behavior of 
referrals and regular clients is the same for both systems. 

First we will discuss the situation where there is no transition between refer¬ 
ral/potential referral and customer/potential customer and thus we set A5 = A7 = 0. 
We consider two cases: the situation where A3 + toA 4 7^ 0 and A 4 + mA 4 7^ 0 (we 
named it “static social network” to reflect the fact that there is no transition be¬ 
tween referrals and regular customers) and the situation where A 4 = A2 = A 4 = 0 , 
corresponding to the case where all potential customers and potential referrals that 
become customers are consequence of referral influence (we named it “word of 
mouth” to emphasise that all marketing efforts are related to referrals). 

We have the following result in the static social network case. 

Theorem 3 (Static social network). The following holds for system 0 with e > 0, 
A2 > 0, A3+mA 4 7^ 0, Ai+mA 4 7^ 0 and A5 = A7 = 0: there is a unique equilibrium 
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solution (R*, C*, Pft, Pq) that is locally asymptotically stable and is given by 


( 7 ) 


R* 


c* 7 _ fl f ary gV «7(g + &) 

2e J £(A 2 +TO fl A 6 )’ 


_ (1 - q:)7(Ai + mM + X 2 R*) 

e{e + /3i + Ai + 771A4 + A2.R*) 


Pr = r !a/e — R* and Pq = 7(1 — ot)/e — C* where 

q _ A3 + TOA4 + /3 2 + £ 
2(A2 + niRXe) 


Define 


( 8 ) 


cry (A2 + m R X 6 ) 
e(e + fa) 


We now consider the word of mouth case. In Figure [2] we show the behavior in 
the plane C — R in both of the regimes of Theorem HJ In Figure [2] we used for the 
plot in the left m = 40, tur = 0, e = 0.01 and A 2 = 10~ 5 , and, for the plot in the 
right, to = 30, tur = 10, e = 0.01 and A 2 = 10 -5 . 


Theorem 4 (Word of mouth). The following statements holds for system (JTJ) with 

£ > 0, A 2 > 0 and Ai = A 3 = A 4 = A 5 = A 7 = 0; 

a) */t ^ 1 then there a unique locally stable equilibrium given by 


{C*,R*,PZ,P*r)= 0,0, 


7(1 — a) 7 a 


b) 


if t > 1 then there are two equilibrium solutions. An unstable equilibrium given 
by 


{ci,r* 1 ,p* C ' 1 ,p*r, 1 )= (0,0, 


7(1 — a) 7a 


and a locally stable equilibrium given by 

( D* P* p* A 

= ( a(l - a)A 2 7 2 (l - 7) 07(1 - 7) (1 - a)(e + ^1)7 £ + / 3 2 A 

Ve 2 (£ + ^1 + A 2 (l - i))’ £ ’ £(£ + /?! + A 2 (l - i))’ A 2 + m R X 6 J 



- \ xWWuWW \ \ \ \ \ \ \ 



Figure 2. Word of mouth: r < 1 and r > 1 
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We next consider a scenario where there is no direct referral influence and thus 
we set A 2 = A 6 = 0. Define 

A 5 [A 7 p + (Ai + 77 iA 4 )<j] A 7 [A 5 g + (A 3 + m\i)p\ 

Kl = '- 7 --WW- - and k 2 = -7-77-7—r-■ 

up(s + (3 2 + A 7 ) vq(e + ft + A 5 ) 

Theorem 5 (No referral influence). For system ([lj with e > 0, A 5 > 0 and 
A 2 = Ag = 0 there is a unique equilibrium solution that is globally asymptotically 
stable and is given by 

(R*,C*,P^,Pc) = (kip, K 2 P, p-Kip, q - H 2 p). 

In Figure [3] we illustrate the behavior of the reduced system in the plane C — R 
for the setting in Theorem [5] In Figure [3] we used A 2 = = 0, m = 40, uir = 0, 

Ai = A 3 = A 4 = A 7 = 0 . 0002 , A 5 = 0.0000018, ft = /3 2 = 0.18, e = 0.01 and, for 
the plot on the left, a = 0, and, for the plot on the right a = 0.5. Note that for 
the figure in the left the number of referrals in the equilibrium point is nonzero, 
although it seems to be (the number of referrals is very low since the there are no 
referrals entering the population). 



Figure 3. Word of mouth: a = 0 and a = 0.5 


3. Simulation 

To obtain a better understanding of the behavior of our model, we assume that, 
for a given corporation, we consider the values for the parameters presented in 
Table [lj and initial conditions Co = 2200, Pc,a = 22000, Rq = 20 and 1ft. 0 = 200. 


Parameter 

Value 

ck, e 

0.01 

7 

eN 0 

Ai,..., A 4 , Ag, A 7 

0.0002 

As 

X 7 R 0 /C 0 

/3i> #2 

0.18 


Table 1. Values of parameters 
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We write No = Co + Pc, o + Ro + Pr,o- The value considered for (3i and fa is 
based on usual assumptions concerning the defection rate |TTJj. The assumptions 
A 5 C 0 = A 7 .R 0 and AsPc.o = MPr,o are made to assure that the underlying social 
network hasn’t an initial tendency to “benefit” any of the four compartments. We 
also consider 7 = No x e so that the total population converges to an equilibrium 
where the total population equals Nq. 

We solved system (JTJ) (named initial ) and system <[ 6 |) (named reduced) with 
MATLAB. In the figures we plot the solution for the initial and the reduced sys¬ 
tems. We considered two sets of values for m and mu, namely (to,?ti_r) = (40,0) 
corresponding to a situation of undifferentiated marketing and ( m,mn) = (30,10) 
corresponding to a situation where some effort is made for attracting referrals. In 
both situations we maintain the same total effort, m + itir = 40, in order to be 
able to compare both cases. 

In Figures|3]and[5]we consider the evolution of customers and referrals in the case 
where marketing is used in an undifferentiated way. We can see that the number 
of customers and referrals decreases in this situation and stabilizes in some lower 
value for both compartments. 
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Figure 4. Regular customers evolution with m = 40 and mu = 0 




Figure 5. Referral evolution with m = 40 and mu = 0 
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In Figures [6] and [7] we now consider the evolution of customers and referrals in 
the case where some marketing effort is used to attract referrals. We can see that 
there is an initial small decrease in the number of customers that is followed rapidly 
by an increase that asymptotically doubles its number. There is also an increase in 
the number of referrals due to the positive value of m#. 


4500 

4000 

3500 

3000 

2500 


-reduced 



(A) t € [0; 7] 



Figure 6. Regular customers evolution with m = 30 and to# = 10 




Figure 7. Referral evolution with to = 30 and to# = 10 

In the Figure [8] we present the evolution of customers and referrals in the case 
where A 2 is reduced to 10 -5 in order to satisfy the condition (|5J) in Theorem ©. 
We can see that, as stated in the theorem, the solutions asymptotically approach 
the same value. In the previous cases, although condition © is not satisfied, 
there is computational evidence that the same happens. Thus we conjecture that 
condition © can be weakened. 


4. Proofs 

4.1. Proof of Lemma [l] Let e > 0. Analysing the direction of the flow on the 
boundary of the set {(x,y, z,w) £ R 4 : x,y,z,w ^ 0} we immediately obtain [aj). 
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(a) Regular customers evolution 



(b) Referrals evolution 


Figure 8 . Asymptotic behavior with m = 30, mu = 10 and A 2 = 
1 x 10“ 5 


Adding the four equations in (HI) and letting N(t) = C(t) + R(t) + Pc it) + Pjj(t), 
we get the differential equation 

(9) TV' = 7 - eTV. 

The general solution of © is N(t) = 7 /e + C e~ £t and thus lim N(t) = 7 /e. In 

£—>•+00 

particular, if N(t) corresponds to an equilibrium solution then N(t) = 7 /e. This 
proves 0 . 


4.2. Proof of Theorem [1} Adding equations for R and Pr and for C and Pc 
in m we get 

J (R + Pr)' = orf — (e + A 7 )(i? + Pr) + A 5 (C + Pc) 

\(C + Pc)' = (1 - +7 - (e + A 5 )(C + P c ) + A 7 (i? + P R ) 

and, setting X = R + Pr and Y = C + Pc, we obtain 

( 10 ) 


X' = orf - (e + A 7 )A + A 5 F 
Y' = (1 - +7 + X 7 X - (e + A 5 )y 


The linear system m has the general solution 


X(t) = 


7 (qg +A 5 ) CiA 7 - C 2 X 5 (£+A5+A7)t (Ci + C 2 ) A 5 t 

e(e + A 5 + A 7 ) A 5 + A 7 A 5 + A 7 


Y (t) = 7((1 - a)s + M) C 2 X 5 - CiA 7 (e+A5+A7)t (C\ + C 2 )X 7 

e(e + A5 + A 7 ) A5 + A 7 A5 + A 7 


and thus 


( 11 ) 


P R (t) = 


7(qg + A 5 ) 

e(e + A 5 + A 7 ) 


0 1 (t) - R(t) 
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with Oi(t) -A- 0 as t -A oo for i = 1,2. Replacing Pr and Pc by the expressions 
in ED in the first two equations in system 0 we obtain 
( 12 ) 

j C' = A7R — (e + / 3 i + A 5 )C + (Ai + U 1 X 4 + A2R) (q + $2(t) — C 1 ) 

[R' = A5C — (£ + /?2 + A7)R + (A3 + U1X4 + (A2 + UIrXo)R) (p + $i(f) — R) 

where p and g are given by ©. 

By Lemma [T] every equilibrium solution must belong to Ao- Thus every equilib¬ 
rium solution must satisfy 


A7R — (£ + /3\ + A 5 )C + (Ai + 171 X 4 + A2R) {q — C) — 0 

A5C — (e + / 3 2 + A7)R + (A3 + ?tiA4 + (A2 + mnXe)R) (p — R ) = 0 


Since A 5 ^ 0 we have 

, , _ (e + /?2 + A 7 )R — (A 3 + TOA 4 + (A 2 + rriRX6)R) (p — R) 

l 14 ; c —-7-• 

Substituting m in the first equation of (fl3l) we obtain equation ©. Each solution 
of © uniquely determines C , R and Pr. Thus, we have at most three equilibrium 
solutions and© follows. 

By (fill for any equilibrium solution (R*,C*,P^, Pq) we have Pr= p — R* and 
P*=q- C*. By UM we obtain © and©. 


4.3. Proof of Theorem [2j Assume that £ > 0 and choose <5 > 0. Additionally, 
define V : (1R + ) 2 —>• R by V(x,y) = \{x 2 + y 2 ) and let (C(t), R(t), Pc(t), PR(t)) 
be a solution of © with initial conditions C(t 0 ) = Co, R(fo) = Ro, Pr{0 ) = Pr, 0 
and Pc{ 0) = Pc, 0 and R a (t)) be a solution of © with initial conditions 

C a {to) = C 0 and R a (to) = Ro- Define x(t) = C(t) - C a {t) and y{t) = R(t ) - R a (t). 
By (fin) and © we have 

x' = C' - C" A 

= X r (R - R 0 ) - (e + Pi + A 5 )(C - C a ) + (Ai + mA 4 )(0 2 (i) - (C - C a )) 

+ qX 2 ^R — R a ) + A 2 02(i)R + X 2 (R a C a — RC) 

(15) = Xry -(e + /?i + A 5 )cc + (Ai + mX 4 )( 0 2 (t) - x) + qX 2 y + X 2 9 2 (t)R 

+ A 2 (Ro(C 0 - C) + (Ro - R)C) 

= —(e -T /3i + A 5 + Ai + 171 X 4 T A 2 Ra)x + (A 7 + qX 2 — CX 2 )y 
+ (A 2 R + (Ai — ?7?.A4))02 (t) 

and 

1/ = R' - K 

= A 5 (C — Co) — (e + @2 + A7 + A3 + 777.A4 + (A2 + rriRXe)p)(R — R a ) 

(16) + (A 2 + itirXg)(R'' — Ro) + ((A 2 + ijirXg)R + A 3 + m.X 4 )di(t) 

= X 5 X - (e + /?2 + A 7 + A 3 + 771 A 4 + (A 2 + rriRXe)p)y 

+ (A 2 + mRXo)(R + R a )y + ((A 2 + itirXo)R + A 3 + m.X4)0i(t) 

By (fill) in the proof of Theorem© there is To ^ 0 such that Rf{(t), R(f) < Rr(£) + 
R(t) < p + 8 and Pc{t),C(t) < Pc{t) + C(t) < q + 5, for all t ^ T 0 . Using this fact, 
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by m and m we have 
j ;V(x,y ) = xx' + yy' 

= — (e + /?i + A5 + Ai + 171X4 + A2 Ra)x^ + (A7 + q\'2 — C\'2)xy 

+ (A2 R + (Ai — m\4))xd2(t) + A5 xy 

— (e + @2 + A7 + A3 + TOA4 + (A2 + niRXo)p)y 2 

+ (A2 + mn\e)(R + R a )y~ + ((A2 + iiirX^R, + A3 + mX^yOiit) 

^ — (e + /?i + A5 + Ai + mX^x^ + (A5 + A7 + 9A2 — CX2)xy 

— (e + P2 + A7 + A3 + 777A4 + (A2 + niRXe)(p + 2 5 ))y 2 + @(i), 

for t ^ T 0 , where 

0(t) — (A 2 -R + (Ai — mX4,))(p + 5)6 2 {t) + ((A 2 + ixirXg)R + A 3 + mXiCj^q + 5)6\{t). 


Thus, using the fact that xy ^ l/2{x 2 + y 2 ), we get 

— P(a;, y) sj — (e + /3i + A 5 + Ai + mXi)x + — (A 5 + A 7 + qX 2 )(x + y ) 


— (e + /?2 + A7 + A3 + mA4 + (A2 + thrXg){jp + 25))y + 0 (t) 

1 1 9 

^ —(s T /Si + ~(As — A 7 ) + Ai + mA 4 — —qX 2 )x 

— (e + /?2 T ~ (A 7 — A 5 ) + A 3 + 1x1X4 + (A 2 + ixirXq)( p + 25) 


- 2 ^A 2 )i/ 2 + 0(t), 

for t ^ To. By J5]), there is M > 0 such that 

d 


dt 


V{x,y)^-MV(x,y) + Q(t). 


Since 0(t) —¥ 0 as t —> 0, there is T\ ^ To > 0 sufficiently large such that, for t ^ Tf, 
we have 0(t) ^ 5. Thus 

d . 


dt 


V(x(t),y(t)) s: -MV(x,y) + 5 


for t > Ti. Therefore, for t^Ti, 


V(x(t),y(t)) < — + P(z(Ti), y (Ti)) e- M(t “ Tl) . 


We conclude that 


lim V(x(t),y(t)) < —. 
t-> + oo M 


Since 5 > 0 is arbitrary, the theorem follows. 

4 . 4 . Proof of Theorem l3l Adding equations for R and Pr and for C and Pc in JT|) we 
get the system 

J (R + Pr)' — ory — e(R + Pr) 

{(C + PcY = (1 - a) 7 - £ {C + Pc) 

and thus 

( R(t) + P R (t) = sa + Cl e~ st 

[C(t) + Pc(t) = £=f^ + C 2 e- st 
Therefore, if (R*, C*, Pr, Pc) is an equilibrium solution, then 
(17) Pr = 07/ e — R* and Pc = (1 — a)7/e — R* ■ 


and 
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Using these expressions and the first two equations in System we obtain 

f—(e + Pi)C* + (Ai + m \ 4 + X 2 R*) -(?*)= 0 

\-(e + p 2 )R* + (A 3 + mA 4 + (A 2 + m R X 6 )R*) - R*) = 0 

and using the second equation in © and dividing by —(A2 + m R Xe) (that is nonzero by 
assumption), we get 


(R *) 2 + 29 R* - 


« 7 (A 3 + m A4) 


= 0 . 


e(A2 + m R Xe) 

Thus we have a unique nonnegative root given by 0. By the first equation in (USD we 
obtain 

(1 — 0)7 Ai + mA4 + A2 R* 


C * = 

by © we conclude that 


e e + Pi -(- Ai T 771A4 -f- A2 R* 


Pr = 


«7 

2 7 


+ 9 - 



c»7(A3 + m\ 4 ) 
e{X 2 + m R \ 6 ) 


and it is easy to check that the expression above is nonnegative. Again by © we get 
Pc = 7(1 — a)/s — C* and it is also immediate that this expression is nonnegative. 

To study the stability of the equilibrium, we consider the Jacobian matrix J at the 
equilibrium. Namely 


J = 


-e - pi A 2 Pc 
0 —£ — A 

Pi -A 2 P* C 

0 A 


Ai T 771A4 H - X 2 R 
0 

—e — Ai — 171X4 — X 2 R* 

0 


0 

B 

0 

—e — B 


where A = p 2 — (A2 + mRXs)P R and B = A3 + mX 4 + (A2 + mRX§)R*. It is easy check 
that the eigenvalues of J are —e, —(A + B + e) and — (Pi + e + Ai + X 4 m — X 2 R*). Thus, 
all eigenvalues have negative real part and we conclude that the equilibrium solution is 
locally asymptotically stable. 


4 . 5 . Proof of Theorem[4j By the second equation in 0 we can see that, if ( R* , C* , Pr, Pq) 
is an equilibrium solution of our system, then 

((A 2 + m R Xe)PR - e - Pi)R* = 0 


and thus R* = 0 or Pr = (e + P)/(X 2 + mRX&). Using the remaining equations and 
condition it is straightforward to obtain the equilibrium points. 

In the equilibrium {R \, Cj, Pr,i,Pc.i), the Jacobian matrix J is given by 

—E — pi A2(l —0)7 /e 0 0 

0 —£ —p 2 + (X 2 +mRXe)'ya/£ 0 0 

Pi —A2(l —0)7 /e — £ 0 

0 P 2 (X 2 + mRXe)-ya/E 0 — e 


We can easily check that the constants —e, —bi — e and 7 q(A2 + XemR^/s — ( p 2 + e) are 
the eigenvalues of J. Therefore if r > 1 the equilibrium is unstable and if r < 1 it is 
locally asymptotically stable. 

In the equilibrium (R 2 , C 2 i Pr, 2i Pc, 2)1 the Jacobian matrix J has the following form 


J = 


-e - Pi 
0 

Pi 

0 


A tPc ,2 X 2 R 2 

0 0 

—X 2 Pq 2 —e — x 2 r 2 

-e 0 


0 

(A2 + mRXe)R 2 

0 

—£ (A2 + mRXe)R 2 


We can easily check that the negative constants — e, —{Pi + e+X 2 R 2 ) and — (X 2 + mRXe)R 2 
are the eigenvalues of J. Therefore the equilibrium is locally asymptotically stable. 
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4 . 6 . Proof of Theorem 03 In our case, equation m in Theorem [T] has a unique solution 
which is nonnegative and is given by 

_ (A3 + m\i)up + (Ai + m.\4)\$q 
uv — A5A7 

Thus, again by Theorem[T| the possible equilibrium solutions are given by (R *, C *, Pr, P£) 
where R* is the constant above, 

_ {s + 02 + A7 + A3 + vn\&)R* — (A3 + mA 4 )p 

° ~ As 

_ uvp(\3 + mA4) + \5vq(\i + mA 4 ) — uvp( A3 + mA 4 ) + A5A7KA3 + mA 4 ) 

X$(uv — A5A7) 

_ (Ai + m\ A )vq + (A 3 + mA 4 )A 7 p 
uv — A5A7 


and 


D » D * up(e + 0 2 + A 7 ) - A 5 (A 7 p+ (Ai + mA 4 )g) 

Pr = V - rt = -T—;- 

UV — A5A7 

vq(e + 0 i + A5) — A7(A5 q + (A3 + mA 4 )p) 

P c = p-C = --- 

uv — A5A7 


Since max{rei, k2} ^ 1 , is immediate that Pr ^ 0 and P ^ 0 . We conclude that there is 
a unique equilibrium solution (R*,C *, Pr, Pc)- 

The Jacobian matrix J at the equilibrium is given by 


J = 


Ai + mA 4 
0 

—e — A5 — Ai — mA 4 
As 


—£ — / 3 \ — A5 A7 Ai + 171X4 0 

A5 —£ — 02 — A7 

Pi 0 

0 02 

We can check that the eigenvalues are all negative and given by — e, — e — A5 — A7 and 


A3 + mA 4 
A7 

—e — A7 — A3 — mA 4 


- a — 2X4171 ± 3/ cr 2 — 4 [(e + 0 i + Ai + As)(e + 02 + A3 + A7) — ArAs]^ . 


Therefore the equilibrium is asymptotically stable. Since the system in this case is linear, 
the equilibrium is globally asymptotically stable. 


5. Conclusions 

We presented and studied in this paper a compartmental model with four compart¬ 
ments to describe the evolution of the number of customers and potential customers of 
some corporation based on the marketing policy of the corporation, determined by the 
effort used in undifferentiated marketing campaigns and in referral directed marketing. 
Apparently the results obtained are reasonable in the sense that the qualitative behavior 
obtained is not going against common sense. The results show that the model works in 
theory in the several scenarios considered, with and without marketing incentives, thus 
relying on normal marketing policies or on incentives to referrals. The model shows that, 
in theory, it is possible to predict the influence referrals can have on their peers based on 
the incentives given to them by the company. The different scenarios also allows to see, in 
a specific period of time, what happens to the number of current and potential customers 
based on the zero incentives policy, and therefore, only based on the natural attractive¬ 
ness power of the referrals (Theorem [ 4 ] with A6 = 0 ); what happens if the company invests 
in marketing but not on incentives to referrals (Theorem O; and also, what happens if 
the company only invests on incentives to referrals (Theorem [I] with A6 7^ 0 ). Based on 
this results, this model allows companies to adjust their marketing policies in order to 
maximize a specific parameter of the model. For instance, the model allows a company 
to estimate the amount of investment necessary to transform an x number of potential 
customers in customers. 
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This study must be followed by work where the model is used in real world problems. 
In fact, comparing the results given by the model and real data would of fundamental 
importance to confirm the usefulness of it. Naturally this is a major task whose feasibility 
will certainly depend on the accuracy in the estimation of parameters. 

We believe that this work opens several possibilities for future studies. For instance, 
it would be interesting to consider versions of this model with time-varying parameters 
to model seasonal phenomena that may occur in some economic activities. It can also 
be of interest to consider age structured populations, to distinguish different consumption 
habits, or to subdivide the universe of referrals, to reflect different aspects that make those 
customers important to the corporation. 
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